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Abstract 

Background: Acute Myeloid Leukaemia (AML) is a highly heterogeneous disease. Studies in adult AML have identified 
epigenetic changes, specifically DNA methylation, associated with leukaemia subtype, age of onset and patient survival 
which highlights this heterogeneity. However, only limited DNA methylation studies have elucidated any associations in 
paediatric AML. 

Methods: We interrogated DNA methylation on a cohort of paediatric AML FAB subtype M5 patients using the lllumina 
HumanMethylation450 (HM450) BeadChip, identifying a number of target genes with p <0.01 and A(3 >0.4 between 
leukaemic and matched remission (n = 20 primary leukaemic, n = 13 matched remission). Amongst those genes 
identified, we interrogate DLEU2 methylation using locus-specific SEQUENOM MassARRAY® EpiTYPER® and an increased 
validation cohort (n = 28 primary leukaemic, n = 14 matched remission, n = 17 additional non-leukaemic and cell lines). 
Following methylation analysis, expression studies were undertaken utilising the same patient samples for singleplex 
TaqMan gene and miRNA assays and relative expression comparisons. 

Results: We identified differential DNA methylation at the DLEU2 locus, encompassing the tumour suppressor microRNA 
miR-15a/16-1 cluster. A number of HM450 probes spanning the DLEU2/Alt1 Transcriptional Start Site showed increased 
levels of methylation in leukaemia (average over all probes >60%) compared to disease-free haematopoietic cells and 
patient remission samples (<24%) (p < 0.001). Interestingly, DLEU2 mRNA down-regulation in leukaemic patients (p < 0.05) 
was independent of the embedded mature miR-1 5a/1 6-1 expression. To assess prognostic significance of DLEU2 DNA 
methylation, we stratified paediatric AML patients by their methylation status. A subset of patients recorded methylation 
values for DLEU2 akin to non-leukaemic specimens, specifically patients with sole trisomy 8 and/or chromosome 
1 1 abnormalities. These patients also showed similar miR-1 5a/1 6-1 expression to non-leukaemic samples, and 
potential improved disease prognosis. 

Conclusions: The DLEU2 locus and embedded miRNA cluster miR-1 5a/l 6-1 is commonly deleted in adult cancers and 
shown to induce leukaemogenesis, however in paediatric AML we found the region to be transcriptionally repressed. 
In combination, our data highlights the utility of interrogating DNA methylation and microRNA in combination with 
underlying genetic status to provide novel insights into AML biology. 
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Background 

Acute myeloid leukaemia (AML) is the third most common 
form of leukaemia in children, typically characterised by the 
rapid proliferation of primitive haematopoietic myeloid 
progenitor cells [1]. Paediatric AML is a highly heteroge- 
neous disease, which presents a major barrier towards 
the development of accurate disease classification, risk 
stratification and targeted therapies within the clinic. 
The French-American-British (FAB) [2] and more re- 
cently World Health Organisation (WHO) [3] classifica- 
tions of leukaemia take into account cell morphology, 
cytogenetic aberrations and common genetic lesions. 
However, not all patients fall into these well-defined cat- 
egories. Additionally, the recurrent chromosomal and 
genetic lesions frequently found in AML fail to induce 
leukaemogenesis and do not explain the recognised 
clinical heterogeneity [4,5]. 

One of the hallmarks of nearly all human cancers is 
the disruption of the epigenetic profile, including gross 
aberrations in DNA methylation. Increasing evidence in 
adult AML has indicated that epigenetic events play 
critical roles in the onset, progression, and outcome of 
AML [6] and may help tailor disease treatment. How- 
ever, the need for similar elucidations in childhood 
disease is paramount. Aberrant methylation of cytosine res- 
idues at palindromic CpG sites (often clustered in dense 
CpG 'islands') near gene promoter regions is widely 
studied in carcinogenesis and haematological malig- 
nancies [6,7]. It is now well established that elevated 
DNA methylation is an important mechanism of gene 
transcriptional inactivation [8,9] and genes such as 
ESR1, IGSF4 and CDKN2B/pl5 are epigenetically si- 
lenced in adult leukaemia [6]. Previous studies have 
subdivided adult AML into 16 epigenetic sub-groups 
based on DNA methylation signatures, correlating with 
patient clinical outcome and distinct from both normal 
haematopoietic cells and normal stages of myeloid 
differentiation [4]. Despite such emerging findings in 
an adult context, the utility of individual DNA methyla- 
tion disruptions in paediatric AML has yet to be fully 
evaluated [6]. 

MicroRNA (miRNA) represent an alternative epigen- 
etic regulator, having been implicated in the regulation 
of critical gene expression networks in plants and ani- 
mals. The role of miRNA in haematopoiesis, cancer and 
disease is also beginning to be appreciated [10,11]. The 
global influence of individual miRNA on the genome is 
difficult to dissect, as miRNA can modulate the expres- 
sion of hundreds of genes, and each gene can harbour 
binding sites for several miRNA [12]. Human miRNA 
are initially transcribed (pri-miRNA), and processed by sev- 
eral complexes to form a 70 bp hairpin-loop (pre-miRNA) 
[13]. After successive enzymatic steps, a miRNA:miRNA* 
complementary duplex is formed where the 'functional' 



strand is combined with RISC (RNA Induced Silencing 
Complex) and Argonaute proteins to guide, and inhibit, 
specific target messenger RNA (mRNA) through base pair 
recognition [14,15]. However, the miRNA transcriptome 
is becoming increasingly complex, emphasised by Next 
Generation Sequencing (NGS) technologies. NGS has 
highlighted that alternate miRNA* transcripts, as well 
as miRNA sequence variants (isomiRs [16]) may play 
a biological role, similar to their canonical miRNA 
relatives [17,18]. 

Links between miRNA deregulation and cancer 
diagnosis were first identified in adult Chronic Lympho- 
blastic Leukaemia (CLL), where the loss or down- 
regulation of tumour-suppressing miRNA cluster miR-lSa/ 
16-1 directly caused leukaemic transformation [19,20]. At 
present, no such association has been identified for child- 
hood leukaemia. The expression of paediatric disease- 
associated miRNA has to date only identified a distinction 
between leukaemia of different lineages and the differenti- 
ation of rearranged AMLs within a limited number of 
cytogenetic subtypes [21,22]. Paediatric MLL can be 
distinguished from others by differentially expressed 
miR-126, miR-146a, miR-181a/b/d, miR-100, miR-21, 
miR-196a/b, miR-29 and miR-125b [21]. However con- 
cordance among studies is often low and the mechan- 
ism of deregulation is often unknown [22,23]. 

Genes encoding miRNA can be regulated epigeneti- 
cally in a similar manner to protein coding genes [22]. 
Studies have demonstrated epigenetically regulated 
miRNA in adult AML, including hypermethylation and 
down-regulation of miR-124a and associated deregula- 
tion of target mRNA EVI1, CEBPA and CDK6 inde- 
pendent of diagnostic cytogenetic subtype (reviewed in 
[22]). Additionally, miR-193a targeting KIT, and miR- 
14b targeting CREB have been identified in adult inves- 
tigations as specifically controlled by DNA methylation 
(reviewed in [22]). However, the identification of DNA 
methylation and miRNA expression connections in 
paediatric leukaemia is lacking. 

Paediatric AML has distinct cytogenetic and clin- 
ical features relative to their adult counterparts 
[5,21,24-26]. Therefore, there is a critical need to im- 
prove our understanding of the biology of childhood 
leukaemia as separate entities, distinct from adult 
disease. Cognisant of this, we aimed to identify differ- 
ential DNA methylation within paediatric AML on a 
genome-scale using defined clinical subtypes and age- 
matched controls. We identified a number of signifi- 
cantly altered DNA methylation loci, with associated 
gene and miRNA expression change, between paediatric 
AML and non-leukaemic counterparts. Specifically we de- 
scribe here the epigenetic deregulation of DLEU2, which 
has associated alterations in downstream miR-lSa/16-1 
miRNA cluster expression. 



Morenos et al. Molecular Cancer 2014, 13:123 
http://www.molecular-cancer.eom/content/1 3/1/1 23 



Page 3 of 1 1 



Results and discussion 

The DLEU2 gene is specifically hypermethylated and 
repressed in paediatric AML subtype M5 

The FAB subtype M5 (monocytic/blastic leukaemia) is a 
distinct subtype with characteristic chromosomal abnor- 
malities including t(8; 16), +8 and various translocations 
involving llq23 and the MLL locus such as t(9;ll), 
t(10;ll), t(ll;19) and others [27]. AML subtype M5 
also has a high proportion of cytogenetically normal 
(CN-AML) patients [27], and those with complex 
karyotypes [28]. A combination of these factors add to 
the overall unfavourable outcome of paediatric M5 
diagnosis [29]. Genome-scale methylation profiling of 
AML M5 bone marrow samples identified 3,352 sig- 
nificantly differentially methylated probes (DMPs) be- 
tween paediatric AML FAB M5 (n = 20) and matching 
non-leukaemic (n = 17) samples. Applying more stringent 
feature selection criteria of an adjusted p-value <0.01 
and A|3 of >0.4 reduced the number of DMPs to 137 
(Additional file 1). 

The list of DMPs included several localising to the long 
non-coding RNA DLEU2 and the embedded miRNA clus- 
ter miR-lSa/16-1 [19,20], previously implicated in adult 



leukaemic [19,20,30]. To date, disruption of this region has 
not been observed in paediatric cancers and as such we 
chose to focus on DLEU2 in subsequent analysis. A total 
of three DLEU2 DMPs had an adjusted p-value <0.01 and 
A|3 of >0.4 (cg05394800, cg20529344, cgl2883980). These 
probes were located within three CpG islands (Chr 13: 50, 
690,000-50,708,000: UCSC human hgl9 assembly) at the 
DLEU2/Altl transcriptional start site (TSS) and north 
shore' (Figure 1), a region up to 2 kb upstream from the 
DLEU2/AM TSS CpG island under investigation [31]. 
Henceforth this will be referred to as the DLEU2 promoter. 

A divergent DNA methylation profile was observed 
between all cases of paediatric AML with elevated 
methylation (mean 64% (48-80% CI) across significant 
probes and non-leukaemic samples (23% methylation; 
7-39% CI of mean) (p < 0.001 for all significant probes) 
(Figure 1; Additional file 2). DNA methylation levels at 
CpG sites interrogated by cg05394800, cg20529344, 
cgl2883980 were confirmed using SEQUENOM EpiTY- 
PER MassArray (Additional file 3 and Additional file 4) in 
the discovery sample set of 20 leukaemic patients and 17 
non-leukaemic controls, as well as a validation set of paedi- 
atric AML M5 patients with heterogeneous cytogenetic 
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Figure 1 Regional interrogation of the DLEU2 gene including significantly differentially methylated probes identified by HM450 
analysis. Bottom: Distribution of HM450 methylation probes across the DLEU2 region of 13q4. Genes located in this region include: DLEUl, 
DLEU2, miR15a/16-l microRNA cluster, TRIM13, KCNRG and miR-3613. Top: HM450 probes identified as significantly differentially methylated 
between paediatric AML (FAB subtype M5) and non-leukaemic specimens have been plotted against genomic location. The leukaemic group 
(n = 16) refers to diagnostic bone marrow from paediatric patients. Non-leukaemic group (n = 1 1) consists of CD sorted cell populations (CD19+, 
CD33+, CD34+, CD45+) and patient remission specimens. Differential methylation is concentrated in the promoter region of the DLEU2/Alt1 long 
transcript, which also falls into the body region of DLEUl within 3 CpG islands (Chr 13: 50, 690,000-50,708,000). Mean methylation (3 values and 
95% CI are shown. Data for leukaemic samples are red and non-leukaemic samples in green. Methylation values range from 0 (0%, no detected 
methylation) to 1.0 (100% fully methylated). Significantly differentially methylated regions between leukaemic and non-leukaemic samples to 
p < 0.001 are highlighted with **. 
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diagnoses (n = 19) (CN-AML (n = 5), MIL (n = 7), RUNX1 
(n = 2), WT1 (n = 1)) (Additional file 5). This analysis con- 
firmed specific DLEU2 promoter hypermethylation in 
paediatric AML patients (p < 0.001) (Additional file 4), with 
additional genes in this region, including DLEU1, hsa-miR- 
lSa/16-1 microRNA cluster, TRIM13, KCNRG and hsa- 
miR-3613 showing no differences in methylation between 
cases and controls (Additional file 6). This highlights the 
potential regional specificity of the observed DLEU2 methy- 
lation change in association with AML. 

We observed a significant down-regulation of DLEU2 
gene expression in paediatric AML (0.07 Fold Change 
(FC); p = 0.014), and a significant inverse correlation be- 
tween promoter DNA methylation and gene expression 
levels (p = 0.0001, Additional file 7). Recent studies in 
adult CLL have also identified a negative correlation be- 
tween DLEU2 promoter methylation (DLEU2/Altl) and 
gene expression [32]. Interestingly there was no change 
in gene expression for any other genes in this region 
(TRIM13 = 1.75 FC; DLEU1 = 1.05 FC. Figure 2A). 

DLEU2 and embedded miR-15a/16-1 are regulated 
independently in paediatric AML 

The miR- 15a/ 16-1 cluster has been described to have 
potent tumour suppressor activity, targeting numerous 
oncogenic and cell cycle regulatory genes [19,33]. The 
cluster is embedded within intron 4 of DLEU2, and has 
been speculated that expression is driven by the DLEU2 
promoter [20,30]. However, we found no correlation be- 
tween DLEU2 expression and miR- 15a/ 16-1 expression 
in paediatric AML, nor down-regulation of the miR-lSa/ 
16-1 miRNA cluster in relation to increasing DLEU2 
promoter DNA methylation (Figure 2B). In contrast to 
previous reports for adult leukaemia, no significant change 
in mature miR- 16 expression was observed between paedi- 
atric AML and control samples, a result we have reported 
elsewhere [34]. These observations are independent of 
the homologous miR- 16 cluster embedded in SMC4 on 
chromosome 3q26, which shows no significant expres- 
sion or DNA methylation changes in association with 
AML (Additional file 8). Similar results have been reported 
for adult CLL [32]. Taken together this data suggests an 
alternate mode of regulation for the miR- 15a/ '16-1 cluster, 
outside of the DLEU2 promoter region. 

Interestingly, recent research has indicated the processing 
mechanisms of miRNA may be affected by cancer, such 
that mature miRNA expression becomes disassociated from 
precursor miRNA levels, and also from the levels of the 
host gene [35]. We found the miR- 15a/ 16-1 primary pre- 
cursor transcript is in fact expressed up to three-fold higher 
in AML patients compared to non-leukaemic counterparts 
(Figure 2B). The increase in primary transcript appears to 
correspond to an increase in mature miR- 16-1* (2.52 FC) 
and miR- 15a* (2.2 FC) expression, and a moderate increase 




Figure 2 Expression analysis of DLEU2 and miR-15/16 in paediatric 

AML Here we interrogate the gene, miRNA and precursor miRNA 

expression in paediatric AML compared to non-leukaemic This interrogation 

includes DLEU2 and embedded miR-15a/16-1 on chromosome 13q4. 

The leukaemic group refers to diagnostic bone marrow from paediatric 

patients. The non-leukaemic group consists of CD sorted cell populations 

(CD19+, CD33+, CD34+, CD45+) and patient remission specimens, and is 

represented by the dashed line at Y= 1. Fold Change (FC) is plotted 

using normalized data and the 2~ AACt method ± SD, and shows the fold 

change calculated from the means of each group. A. Gene expression 

including DLEU1, DLEU2 and TRIM13 in paediatric AML (n = 10) compared 

to non-leukaemic (n = 13) expression. DLEU2 shows a significant 

down-regulation in AML compared to non-leukaemic expression 

(0.07 FC, p = 0.014 represented by **), however there is no significant 

change in expression for TRIM13 or DLEU1. B. Mature microRNA 

expression, including primary precursor transcript (PRI) and alternate 

miRNA expression (*), from the miR-15a/16-1 miRNA cluster embedded 

within DLEU2 for paediatric AML (n = 28, including the 1 0 used in 

Figure 2A) compared to non-leukaemic specimens (n = 30, including 

the 13 used in Figure 2A). Here the miR-1 5a/1 6-1 PRI transcript is 

3.03-fold higher in expression compared to non-leukaemic expression. 

Additionally, miR-1 6-1* (2.52 FC), miR-1 5a* (2.24 FC) and miR-1 5a (1.5 FC) 

also show increases in expression in paediatric AML. No significant 

change in expression was observed for miR-1 6 (0.94 FC) in paediatric 

AML compared to non-leukaemic expression. 
^ J 
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in miR-15a (1.5 FC). However, individual patient samples 
show differential degrees of over-expression of one, or 
all, of the mature species despite the common precursor 
(Additional file 9). 

Canonical miR-15a/16-l microRNA species target nu- 
merous oncogenic and integral cell cycle regulatory 
genes such as BCL-2, MCL-1, CCND1, CDK6, BMI-1, 
RASSF5, IGSF4, c-MYB and WNT3A [33,36-42]. To iden- 
tify the potential significance of increasing alternate 
miR- 15a/ 16-1 transcripts over their canonical counter- 
parts in paediatric AML, we investigated miRNA target 
genes using prediction tools, and undertook gene ontol- 
ogy analysis. We found miR-16-1* is three-fold enriched 
for targeting components of the RNA processing and 
splicing machinery such as hnRNP regulators, SRRM1 
and NUDT2. Down-regulation of these genes by increased 
expression of a targeting miRNA may be contributing to 
miRNA/host gene disassociation in DLEU2/miRNA-15a/ 
16-1 expression (Additional file 10). One target of miR-15a* 
was found to be ADORA2A, an important G-protein recep- 
tor critical in tissue-specific and systemic inflammatory 
responses [43,44]. We have identified the down-regulation 
of ADORA2A expression within paediatric AML, not me- 
diated through DNA methylation, but potentially through 
the up-regulation of the miR-15a/16-l cluster (data not 
shown). Moreover, miR-15a* and miR-16-1* are 60-90 fold 
enriched for targeting genes involved in the intrinsic 
apoptotic pathway such as BCL-2 Lll, PPIF and DNM1L 
(Additional file 10). The inhibition of these genes to act 
upon the cell, mediated through targeting miRNA, can po- 
tentially encourage continued cell growth and perpetuation 
of AML, through negating apoptotic signalling. The over 
expression of alternate mature miRNA in paediatric AML 
may also have consequences for integral cell cycle path- 
ways, but through different mechanisms to their canonical 
counterparts. 

DLEU2 interrogation identifies a novel subclass of 
paediatric AML 

In addition to assessing the association of DLEU2 promoter 
DNA methylation and miR- 15a/ 16-1 miRNA expression 
with paediatric AML, we also assessed correlations with 
distinct clinical and diagnostic variables. We found no sig- 
nificant association between DLEU2 methylation/expres- 
sion with sex, age of disease onset or relapse status nor 
common gene abnormalities such as FLT3 or MLL, as has 
been found in a number of other association studies [6,45]. 
DLEU2 methylation clustering based on diagnostic FAB 
subtype revealed a wide range of values (M5a: 16-91%; 
M5b: 48-90%; 'Other' subtypes: 51-93%. Additional file 11), 
indicating mean methylation in isolation does not stratify 
paediatric AML subtype. Additionally, we found there to be 
no correlation between percentage of leukaemic blasts at 



patient diagnosis and the DNA methylation status at the 
DLEU2 promoter region (data not shown). 

Previous adult leukaemia studies have identified novel 
methylation subgroups, independent of cytogenetic pro- 
file [4,46]. We identified a subset of patients displaying 
low DLEU2 promoter methylation comprising patients 
with trisomy 8 and/or chromosome llq abnormalities, 
independent of MLL rearrangements (separated into a t 
(11)/ +8 subgroup - Figure 3 A, Additional file 11). This 
delineation identified t(ll)/+8 subgroup as having com- 
parable DLEU2 methylation levels to non-leukaemic 
specimens (Figure 3B), and significantly lower methyla- 
tion compared to the remaining paediatric AML samples 
(M5a/b p < 0.05; M1/M2/M4 p < 0.0001 Figure 3B). Classi- 
fying t(ll)/+8 patients as a separate subgroup further 
delineated M5a, M5b and M1/M2/M4 AML patients from 
non-leukaemic counterparts according to average methyla- 
tion status (M5a MM = 0.81, M5b MM = 0.78, M1/M2/M4 
MM = 0.87. p < 0.0001. Figure 3B). 

Cases carrying llq rearrangements are the most het- 
erogeneous of paediatric AML [27,29,47], linked to 
50-104 translocation fusion partners to date [29,46,48]. 
Trisomy 8 is a frequently reported aberration in adult 
and paediatric AML [29]. However little is known about 
the gain of chromosome 8 in isolation and its relation- 
ship to disease onset. As such this has been speculated 
to be a disease modulating secondary event [49]. Lower 
DNA methylation in the t(ll)/+8 subgroup may poten- 
tially confer a better prognosis, as it is well documented 
that hypermethylation of the DLEU2 region is involved 
in leukaemic transformation in adults [19,20,30]. Ana- 
lysis of the t(ll)/+8 subgroup in isolation also revealed 
no significant miRNA expression differences relative to 
non-leukaemic samples (Figure 4A), and it is well docu- 
mented that miR- 15a/ 16-1 abnormalities are also associ- 
ated with leukaemic transformation in adults [19,20]. 
We additionally identified a trend towards decreased risk 
of relapse in t(ll)/+8 subgroup compared to other sub- 
types, and a trend towards better survival outcomes 
(Additional file 12). These analyses combined may elu- 
cidate a connection between DLEU2 promoter DNA 
methylation (Figure 3B), miRNA expression (Figure 4A) 
and prognostic outcomes for this subgroup of paediat- 
ric AML (Additional file 12). 

The classification of t(ll)/+8 cases as an independent 
subgroup revealed a significant increase in pri-miR-15a 
expression in traditional FAB subtype M5b from non- 
leukaemic (20.29-fold, p< 0.001; Figure 4C) and from 
M5a/Ml/M2/M4 (Figure 4B). The increase translates to 
a >2-fold increase in miR- 15a, 4.87-fold increase in miR- 
15a*, and 9.86-fold increase in miR-16-1* expression. 
Based on DLEU2 DNA methylation status, in conjunc- 
tion with miRNA expression analysis, we speculate that 
pri-miR-15a expression alone may be a useful biomarker 
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Figure 3 Interrogation of paediatric AML by clinically defined cytogenetic and FAB subtypes alongside DLEU2 promoter methylation. 

Paediatric AML diagnostic cytogenetic and subtyping analyses were specifically investigated, including those with known gene abnormalities 
and CN-AML cases. We investigate here the DNA Methylation of DLEU2/Altl promoter region to 95% CI. DNA methylation values range from 0.0 
(0% no detected methylation) to 1.0 (100% fully methylated). The leukaemic group refers to diagnostic bone marrow from paediatric patients. 
Non-leukaemic group consists of CD sorted cell populations (CD19+, CD33+, CD34+, CD45+) and patient remission specimens. A. DNA methylation for 
DLEU2 HM450 probes (eg 12883980, cg20529344, cg5394800) according to cytogenetic type, showing heterogeneous DNA methylation outcomes. Of note, 
a subset of patients with observable chromosome 1 1 and trisomy 8 abnormalities have a reduced DNA methylation compared to all other 
AML abnormalities. B. A subset of DLEU2 DNA methylation results for paediatric AML chromosome 1 1 and trisomy 8 abnormality patients 
from Figure 3A fall into a range akin to non-leukaemic specimens. Grouping these patients, regardless of clinical subtyping (M5a, M5b or 
M1/M2/M4), forms an additional sub-group, denoted here as 't(1 1)/+8'. The DNA methylation values obtained for probe eg 1 2883980 are 
depicted, as an indication of the methylation at the DLEU2 promoter. t(1 1)/+8 are not differentially methylated compared to non-leukaemic 
specimens (t(1 1)/+8 sub-group Mean Methylation (MM) = 0.38, non-leukaemic MM = 0.32. p = 0.22). M5a, M5b and M1/M2/M4 AMLs (with the removal of 
t(1 1)/+8 patients) have become increasingly hypermethylated compared to non-leukaemic methylation (M5a MM = 0.81, M5b MM = 0.78, M1/M2/M4 
MM = 0.87. p < 0.0001), and significantly different from t(1 1)/+8 (M5a/b p < 0.05; M1/M2/M4 p < 0.0001). 



to distinguish M5b FAB subtype from all other AML 
subtypes. 

Conclusions 

Previous research within paediatric AML has shown that 
well-defined cytogenetic subgroupings exhibit a wide 
range of genomic and epigenomic heterogeneity. Linking 
specific epigenetic features to clinical parameters has the 
potential to identify pathological drivers of disease and 
develop enhanced molecular approaches for diagnosis, 
prognosis and refinement of treatment. Profiling epigenetic 
regulators has provided clinically relevant biomarkers for 
adult cancers; however generally the same cannot be said of 
childhood cancers. Our analysis has identified hypermethy- 
lation induced down-regulation of the DLEU2 gene in 
paediatric AML. The related expression changes of the 
embedded miR-lSa/16-1 microRNA cluster in paediatric 
AML has the potential to contribute to leukaemic trans- 
formation, with a switch from canonical to alternate 
miRNA family members, which in turn may modulate the 
expression of downstream regulatory genes. Treating 
paediatric AML model systems with epigenetic modifying 
drugs will allow a more comprehensive analysis towards 
the clinical applications and potential patient-focussed 
therapeutic interventions for children. Our results high- 
light the need for further specific interrogation of paediatric 
AML subtypes as distinctive biological entities, separate 



from adult disease. Further studies utilising larger patient 
cohorts are required to explore the complex interplay 
between the epigenetic regulation of genes harbouring 
microRNA, taking into account alternate miRNA tran- 
script expression. 

Methods 

Samples 

This study was approved by the Royal Children's Hospital 
(RCH), Melbourne, Ethics Committee (HREC reference 
#27138E). Samples used consisted of snap frozen bone mar- 
row specimens and archived bone marrow biopsies taken 
at diagnosis, patient remission/follow-up or at relapse from 
paediatric acute myeloid leukaemia (AML) cases. The diag- 
nosis of AML was established according to the criteria of 
the French- American-British (FAB) classification by stand- 
ard morphological and cytological methods. The evaluation 
for mutations in MLL, FLT3, RUNX1 and WT1 genes were 
assessed in a small number of patients as a part of the ini- 
tial assessment. No patient recorded any chromosome 13 
deletions or abnormalities listed in their clinical diagnoses. 
The median percentage of leukaemic blasts at patient diag- 
nosis is 88% (67-93 95% CI). 

Patient samples used in our study consisted of archived, 
air-dried bone marrow smear slides. The utility of these 
samples for DNA methylation and miRNA expression ana- 
lysis has been outlined previously [34,50]. Cryogenically 
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Figure 4 Gene, primary, miRNA and miRNA* expression for paediatric AML defined through clinical classification and DLEU2 
Methylation subtyping. Mature microRNA expression, primary precursor transcript (PRI) and alternate miRNA isoform expression (*) from the 
miR-!5a/l6-1 miRNA cluster embedded within DLEU2 for paediatric AML patients (n = 26: 12 M5a, 5 M5b, 4 M1/M2/M4, 5 t(1 1)/+8 sub-group) all 
compared to non-leukaemic specimens (n = 30). Leukaemic groups refer to diagnostic bone marrow from paediatric patients. Non-leukaemic 
group consists of CD sorted cell populations (CD19+, CD33+, CD34+, CD45+)and patient remission specimens. Linear Fold Change (FC) is plotted 
using normalized data and the 2" AACt method ± SD, and shows the fold change calculated from the means of each group. DLEU2 gene expression 
is down-regulated in all subtypes. A. Fold change in expression comparing non-leukaemic to t(1 1)/+8 subtype. No significant differences in RNA 
expression are observed between non-leukaemic specimens and this subgroup. B. Fold change in expression comparing patients from subtype 
M5a to M1/M2/M4. DLEU2 is down-regulated in all subtypes (as previously described in Figure 2), and primary precursor for miR-15a also appears 
down-regulated (non-significant). M1/M2/M4 groupings do not show any mature miRNA expression changes (defined as > 2-fold difference from 
non-leukaemic). Subtype M5a shows a 2.1 -fold (±0.8 SD) increase in miR-15a* and a 2.81 -fold (±0.4 SD) increase in miR-16-1* C. FAB subtype M5b 
shows up-regulation of miR-15a PRI compared to non-leukaemic expression (20.29-fold (±0.6 SD) p < 0.001), with miR-15a* (4.87-fold ±1 SD) and 
miR-16-1* (9.86-fold ±1 SD) also up-regulated. M5b additionally shows a significant up-regulation of miR-15a PRI compared to t(1 1)/+8 sub-group 
samples (p < 0.05), and also from M5a and M1/M2/M4 sub-groups (p < 0.001). 



frozen patient bone marrows were also used where 
available. All patients were <18 years of age. We chose 
to focus specifically on the FAB subtype M5 (a and b) 
as it is a common subtype. Our cohort included com- 
parable ratios between males and females, and also 
similar numbers of children that relapsed or not. 
(Additional file 5). Control samples from bone marrow 
of unrelated and unaffected children were analysed in 
parallel, as well as multiple cell lines, including adult 
leukaemia (REH [51], CCRF-CEM [52]), paediatric 
AML (Kasumi-1 [53], THP-1 [54], MV-4-11 and AML- 
193 [55]) obtained from American Type Culture Col- 
lection (ATCC) and subjected to characterisation using 



the ATCC Proficiency Standard program. Fluorescent 
Activated Cell Sorted (FACS) isolated haematopoietic 
progenitor cell populations from unrelated and unaffected 
paediatric donors (CD19+, CD33+, CD34+ and CD45+ 
populations, from herein known as CD sorted cells), 
were also used in this study. Our 'non-leukaemic' 
group consisted of all individual CD sorted cell popula- 
tions as well as non-leukaemic remission and follow- 
up slides and bone marrow from paediatric patients. 
We compared the DLEU2/AIU DNA methylation for 
all 'non-leukaemic' specimens, and found no signifi- 
cant differences across these non-leukaemic samples 
(Additional file 13). 
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DNA extraction, quality control and methylation analysis 

Genomic DNA was extracted using the phenol/chloro- 
form method and bisulphite converted in accordance 
with manufacturers protocols, as reported previously [23]. 
DNA extracts were checked for quality and quantity using 
a NanoDrop® ND-1000 spectrophotometer (Thermo Fisher 
Scientific Inc., Scoresby, Victoria, Australia). All DNA was 
stored at -80°C. 

Bisulphite conversion of genomic DNA was performed 
using the MethylEasy Xceed Bisulfite Modification Kit 
(Human Genetic Signatures, Sydney, AUST). The con- 
verted samples were processed by the Australian Genome 
Research Facility (AGRF, Melbourne, Australia) and ana- 
lysed using the Illumina HumanMethylation450 (HM450) 
Bead Chip arrays according to manufacturers protocol. 
Illumina Genome Studio software was used to extract the 
raw M-values and probe intensities for downstream pro- 
cessing. Samples with a probe p-detection value of p <0.05 
were retained. Data were quantile normalized using lumi 
[56] and analysed using LIMMA [57]. All analysis was 
performed using the R statistical software package on 
M-values that were converted to (3-values (0 = unmethy- 
lated, 1 = fully methylated) for reporting and biological 
interpretation. To eliminate sex bias, probes hybridizing to 
sex chromosomes were filtered out, leaving 366,553 probes 
common to all samples in the final dataset. Further filtering 
was based on the degree of difference in |3 between sample 
groups, indicated as Ap calculated using the formula: 

Afi = meany5value(Leukaemic)-meanySvalue(Non-Leukaemic) 

SEQUENOM MassARRAY® EpiTYPER® was used to 
measure locus-specific methylation. Sixty-one samples 
were analysed, including the 42 analysed on the HM450 
platform (Additional file 5), an additional 3 cell lines, 
8 leukaemic and 8 non-leukaemic whole bone marrow 
aspirate samples. Primers for analysis were designed using 
SEQUENOM EpiDesigner software (www.epidesigner.com) 
and the sequences for these are listed in Additional file 14. 
Gene ontology and pathway analysis of genes associated 
with significantly altered DNA methylation probes were 
analysed through the use of IPA (Ingenuity® Systems, 
www.ingenuity.com) and GOrilla [58]. 

RNA extraction, quality control and expression analysis 

Small RNA extraction was performed using TRIzol® 
(Ambion®) for patient bone marrow samples, mononuclear 
cells and cell lines in accordance with the manufacturers 
instructions, or the High Pure miRNA Isolation Kit 
(Roche) for archived slide samples as previously described 
[34]. Before RNA extraction of fresh bone marrow aspi- 
rates, samples were processed using Ficoll-Paque™ (GE 
Healthcare, Piscataway USA) to isolate the mononuclear 
cell population. This was immediately cryo-frozen or stored 



in RNAlatef (Ambion® by Life Technologies, Mulgrave, 
Victoria Australia) for later extraction. The concentration 
and purity of all RNA samples was assessed using the 
NanoDrop® ND-1000 spectrophotometer (Thermo Fisher 
Scientific Inc., Scorsby, Victoria, Australia). All RNA was 
stored at -80°C. 

Seventy samples were used in the interrogation of 
miRNA expression, 64 of which were also used for HM450 
and SEQUENOM methylation analysis (Additional file 5). 
TaqMan® MicroRNA Reverse Transcription kit and single- 
plex TaqMan @ microRNA Assays (Applied Biosystems, Life 
Technologies) (assays listed in Additional file 15) were 
utilised according to the manufacturer s instructions 
before routine quantitative real-time PCR (qRT-PCR) 
was performed using Applied Biosystems 7300 Se- 
quence Detection System. A subset of 29 high quality 
samples which included all cell lines as previously de- 
scribed, CD19+, CD33+, CD34+, CD45+ as well as 8 
leukaemic and 9 non-leukaemic patient whole bone mar- 
rows, were used for gene and pri-miR expression analysis. 
This was due to longer RNA species often being degraded 
in archived specimens [59,60]. The Superscript® VILO™ 
cDNA Synthesis Kit (Life Technologies) and TaqMan® gene 
expression assays were used for gene/pri-miR expression as 
per manufacturers instructions (Additional file 15 and 
Additional file 16) before performing qRT-PCR analysis ac- 
cording to the manufacturers instructions. All qRT-PCR 
samples were analysed in duplicate. 

Expression data was normalized using the synthetic small 
nucleolar RNA, C/D Box 44 (RNU44) and glyceraldehyde- 
3-phosphate dehydrogenase (GAPDH) for gene/pri-miR 
expression; and RNU44 and hsa-miR-26b for miRNA ana- 
lysis, previously identified as having the greatest stability in 
our samples [34]. Data were analysed using DataAssist™ 
software (v.3.0, Applied Biosystems, Life Technologies: 
www.lifetechnologies.com/au/en/home/technical-resources/ 
software-downloads/dataassist-software), and all p-values 
were adjusted using Benjamini-Hochberg False Discov- 
ery Rate. 

All analyses used the combined average of non- 
leukaemic primary bone marrow tissue and CD sorted 
non-leukaemic cell expression values as the Reference 
Group (calibrator). Fold Change (FC) was calculated 
using the Livak method of 2t AAC [61] plotting fold 
change and associated p-values. The fold change re- 
ported here is the difference of the means of each 
group, and a fold change of >2 between disease and 
non-leukaemic groups were considered noteworthy. 
microRNA gene target prediction was assessed using 
microRNA.org [62,63] and miRWalk [64]. The top 500 
target genes for each miRNA (Additional file 17) were 
used for gene ontology/pathway analysis (Additional 
file 10) by IPA (Ingenuity® Systems, www.ingenuity.com) 
and GOrilla [58]. 
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Statistical methods 

Where groups were compared with more than one vari- 
able measured (non-repeated measures), a Two-Way 
ANOVA with Bonferroni Post-testing was performed. 
Where groups were compared with only one variable 
measured, a Kruskal-Wallis test was performed with 
Dunns Post-testing. These analyses were undertaken on 
normalized data. Correlation assessments utilised Spear- 
man Rank tests. Two-tailed p-values were utilised for 
assessment of statistical significance, and results with 
p-values p < 0.05 were considered to be statistically sig- 
nificant. Survival and relapse curves were estimated by 
the Kaplan-Meier method. Statistical analyses were per- 
formed with GraphPad PRISM (Version 5) for Windows, 
(GraphPad Software, USA, www.graphpad.com). CI re- 
fers to 95% Confidence Intervals. 

Additional files 



Additional file 17: The top 500 potential target genes identified for 
each differentially expressed miR-15a/16-l cluster microRNA 
member in paediatric AML. 



Abbreviations 

AML: Acute myeloid leukaemia; CN-AML: Cytogenetically normal acute 
myeloid leukaemia; DLEU2: Deleted in lymphocytic leukaemia 2; 
DMP: Differentially methylated probe; FAB: French-American-British 
leukaemia classification method; HM450: lllumina humanmethylation450 
bead chip array; miR: microRNA/miRNA; Pri-miRNA: Primary microRNA; 
TSS: Transcriptional start site. 

Competing interests 

The authors declare that they have no competing interests. 
Authors' contributions 

MH and MP-B carried out cell sorting, extractions and the accumulation of 
samples. JN assisted in accumulation of samples, extractions, conversions 
and assay running, as well as preliminary data analysis. FM and NG facilitated 
sample procurement and conceived the study. ZC participated in statistical 
and data analysis, study concepts and helped to draft the manuscript. RS 
and NW conceived the study, acquired the samples, interpreted the data, 
participated in drafting the manuscript and provided critical revisions for 
the approved final version. LM designed the study, acquired the samples, 
organised sample processing, participated in sample processing, undertook 
all data analysis and statistical evaluations and drafted the manuscript. All 
authors read and approved the final manuscript. 

Acknowledgements 

We thank Dr Elizabeth Algar and Priscilla Siswara for access to flash frozen 
bone marrow and gene abnormality identification, and Dr. Pei Tian for 
facilitating sample procurement. Dr. Jovana Maksimovic we thank for HM450 
data analysis. 

Financial support 

This work was supported by the Victorian Cancer Agency Grant and National 
Health & Medical Research Council (NHMRC), Australia [Project Grant number 
607306] (N.W.); MCRI is supported by the Victorian Government's Operational 
Infrastructure Support Program; the National Health & Medical Research 
Council Dora Lush Postgraduate Scholarship to LM.; NHMRC Senior Research 
Fellowship to R.S.;The Leukaemia Foundation of Australia Phillip Desbrow 
Postdoctoral Fellowship (N.W.) and Project Grant in Aid( N.W. & R.S.); N.W. 
was also supported by My Room and the Children's Cancer Centre 
Foundation. 

Author details 

dancer & Disease Epigenetics, Murdoch Childrens Research Institute, 
Melbourne, Victoria, Australia, department of Paediatrics, The University of 
Melbourne, Melbourne, Victoria, Australia. 3 Cancer Research, Murdoch 
Childrens Research Institute, Melbourne, Victoria, Australia. 4 BMDI Cord Blood 
Bank, Murdoch Childrens Research Institute, Melbourne, Victoria, Australia. 
5 Current address: Ludwig Institute of Cancer Research, Olivia Newton John 
Cancer and Wellness Centre, Austin Hospital, Heidelberg, Victoria, Australia. 

Received: 4 December 2013 Accepted: 15 May 2014 
Published: 24 May 2014 

References 

1. Faderl S, Talpaz M, Estrov Z, O'Brien S, Kurzrock R, Kantarjian HM: The 
biology of chronic myeloid leukemia. N Engl J Med 1 999, 341 :1 64-1 72. 

2. Bennett JM, Catovsky D, Daniel M-T, Flandrin G, Galton DAG, Gralnick HR, 
Sultan C: Proposals for the classification of the acute leukaemias. Br J 
Haematol 1976, 33:451-458. 

3. Vardiman JW, Thiele J, Arber DA, Brunning RD, Borowitz MJ, Porwit A, Harris 
NL, Le Beau MM, Hellstrom-Lindberg E, Tefferi A, Bloomfield CD: The 2008 
revision of the World Health Organization (WHO) classification of 
myeloid neoplasms and acute leukemia: rationale and important 
changes. Blood 2009, 1 14:937-951. 

4. Figueroa ME, Lugthart S, Li Y, Erpelinck-Verschueren C, Deng X, Christos PJ, 
Schifano E, Booth J, van Putten W, Skrabanek L, Campagne F, Mazumdar M, 



Additional file 1: HM450 methylation heat map of the top 137 
significantly differentially methylated probes between paediatric 
AML (FAB subtype M5) compared to non-leukaemic controls. 

Additional file 2: Paediatric AML DNA methylation results compared 
to non-leukaemic methylation for the differentially methylated probes 
identified using HM450. 

Additional file 3: DNA methylation cross-platform validation between 
HM450 and SEQUENOM MassARRAY® EpiTYPER®. 

Additional file 4: CpG locations captured by SEQUENOM MassARRAY® 
EpiTYPER®, including the HM450 probes. 

Additional file 5: Patient information for all paediatric AML patients 
used for HM450 DNA methylation interrogation, SEQUENOM 
MassARRAY® EpiTYPER® validation and gene expression analysis. 

Additional file 6: HM450 probe methylation values for the 13q4 
region for paediatric AML (FAB M5) compared to non-leukaemic 
controls. 

Additional file 7: Correlation between DLEU2/Altl promoter DNA 
methylation and DLEU2 gene expression. 

Additional file 8: HM450 probe methylation and expression 
dynamics of paediatric AML compared to non-leukaemic controls 
for the alternate miR-16 cluster, miR-15b/16-2 embedded in SMC4 
on chromosome 3. 

Additional file 9: Comparison of individual patient miRNA miR-15a/ 
16-1 cluster expression. 

Additional file 10: Top Gene Ontology terms using GOrilla gene 
ontology identification tool. 

Additional file 11: DLEU2 promoter DNA methylation divided into 
paediatric AML patient clinical diagnostic observations. 

Additional file 12: Survival and relapse curves of paediatric AML, 
including a newly defined sub-group based on DLEU2 promoter 
DNA methylation. 

Additional file 13: DLEU2 DNA methylation comparison of the three 
'non-leukaemic' sample types. 

Additional file 14: Probes used in SEQUENOM MassARRAY® 
EpiTYPER® validation analysis of DLEU2/Altl HM450 promoter 
probes located on Chromosome 13. 

Additional file 15: miR-15a/16-1 cluster TaqMan® microRNA and 
primary precursor expression assays utilized, including 
Normalization references. 

Additional file 16: TaqMan® gene expression assays utilized. 



Morenos et al. Molecular Cancer 2014, 13:123 
http://www.molecular-cancer.eom/content/1 3/1/1 23 



Page 10 of 11 



Greally JM, Valk PJM, Lowenberg B, Delwel R, Melnick AM: DNA methylation 
signatures identify biologically distinct subtypes in acute myeloid 
leukemia. Cancer Cell 2010, 17:13-27. 

5. Pui C-H, Carroll WL, Meshinchi S, Arceci RJ: Biology, risk stratification, 
and therapy of pediatric acute leukemias: an update. J Clin Oncol 201 1, 
29:551-565. 

6. Alvarez S, Suela J, Valencia A, Fernandez A, Wunderlich M, Agirre X, Prosper 
F, Martin-Subero Jl, Maiques A, Acquadro F, Rodriguez Perales S, Calasanz 
MJ, Roman-Gomez J, Siebert R, Mulloy JC, Cervera J, Sanz MA, Esteller M, 
Cigudosa JC: DNA methylation profiles and their relationship with cytogenetic 
status in adult acute myeloid leukemia. PLoS One 2010, 5:e12197. 

7. Gardiner-Garden M, Frommer M: CpG islands in vertebrate genomes. 
JMol Biol 1987, 196:261-282. 

8. Galm 0, Wilop S, Luders C, Jost E, Gehbauer G, Herman JG, Osieka R: 
Clinical implications of aberrant DNA methylation patterns in acute 
myelogenous leukemia. Ann Hematol 2005, 84:39-46. 

9. Bird AP: DNA methylation patterns and epigenetic memory. Genes Dev 
2002, 16:6-21. 

10. Bissels U, Bosio A, Wagner W: MicroRNAs are shaping the hematopoietic 
landscape. Hoemotologico 2012, 97:160-167. 

11. Calin GA, Croce CM: MicroRNA signatures in human cancers. Nat Rev 
Cancer 2006, 6:857-866. 

12. Krek A, Grun D, Poy MN, Wolf R, Rosenberg L, Epstein EJ, MacMenamin P, 
da Piedade I, Gunsalus KC, Stoffel M, Rajewsky N: Comninatorial microRNA 
target predictions. Nat Genet 2005, 37:495-500. 

13. Calin GA: MicroRNAs and cancer: what we know and what we still have 
to learn. Genome Med 2009, 1:78. 

14. Tsai LM, Yu D: MicroRNAs in common diseases and potential therapeutic 
applications. Clin Exp Pharmacol Physiol 2010, 37:102-107. 

15. Zhao G, Yu D, Weiss MJ: MicroRNAs in erythropoiesis. Curr Opin Hematol 
2010, 17:155-162. 

16. Marti E, Pantano L, Banez-Coronel M, Llorens F, Minones-Moyano E, Porta S, 
Sumoy L, Ferrer I, Estivill X: A myriad of miRNA variants in control and 
Huntington's disease brain regions detected by massively parallel 
sequencing. Nucleic Acids Res 2010, 38:7219-7235. 

17. Llorens F, Bahez-Coronel M, Pantano L, Antonio del Rio J, Ferrer I, Estivill X, 
Marti E: A highly expressed miR-101 isomiR is a functional silencing small 
RNA. BMC Genomics 2013, 14:104. 

18. Neilsen CT, Goodall GJ, Bracken CP: IsomiRs - the overlooked repertoire in 
the dynamic microRNAome. Trends Genet 2012, 28:544-549. 

19. Calin GA, Dumitru CD, Shimizu M, Bichi R, Zupo S, Noch E, Alder H, Rattan S, 
Keating MJ, Rai K, Rassenti L, Kipps TJ, Negrini M, Bullrich F, Croce CM: 
Frequent deletions and down-regulation of micro-RNA genes miR15 and 
miR16 at 13q14 in chronic lymphocytic leukemia. Proc Natl Acad Sci USA 
2002, 99:15524-15529. 

20. Klein U, Lia M, Crespo M, Siegel R, Shen Q, Mo T, Ambesi-lmpiombato A, 
Califano A, Migliazza A, Bhagat G, Dalla-Favera R: The DLEU2/miR-15a/l6-l 
cluster controls B cell proliferation and its deletion leads to chronic 
lymphocytic leukemia. Cancer Cell 2010, 17:28-40. 

21 . Daschkey S, Rottgers S, Giri A, Bradtke J, Teigler-Schlegel A, Meister G, Borkhardt 
A, Landgraf P: MicroRNAs distinguish cytogenetic subgroups in pediatric AML 
and contribute to complex regulatory networks in AML-relevant pathways. 
PLoSOne 2013, 8:e56334. 

22. Agirre X, Martinez-Climent JA, Odero MD, Prosper F: Epigenetic regulation 
of miRNA genes in acute leukemia. Leukemia 2012, 26:395-403. 

23. Chatterton Z, Morenos L, Saffery R, Craig JM, Ashley D, Wong NC: DNA 
methylation and miRNA expression profiling in childhood B cell acute 
lymphoblastic leukemia. Epigenomics 2010, 2:697-708. 

24. Meshinchi S, Alonzo TA, Stirewalt DL, Zwaan M, Zimmerman M, Reinhardt D, 
Kaspers GJL, Heerema NA, Gerbing R, Lange BJ, Radich JP: Clinical implications 
of FLT3 mutations in pediatric AML. Blood 2006, 108:3654-3661. 

25. Juhl-Christensen C, Ommen HB, Aggerholm A, Lausen B, Kjeldsen E, Hasle 
H, Hokland P: Genetic and epigenetic similarities and differences 
between childhood and adult AML. Pediatr Blood Cancer 2012, 58:525-531. 

26. Liang D-C, Liu H-C, Yang C-P, Jaing T-H, Hung l-J, Yeh T-C, Chen S-H, Hou 
J-Y, Huang Y-J, Shih Y-S, Huang Y-H, Lin T-H, Shih L-Y: Cooperating gene 
mutations in childhood acute myeloid leukemia with special reference 
on mutations of ASXL1, TET2, IDH1, IDH2, and DNMT3A. Blood 2013, 
121:2988-2995. 

27. Tallman MS, Kim HT, Paietta E, Bennett JM, Dewald G, Cassileth PA, 
Wiernik PH, Rowe JM: Acute Monocytic Leukemia (French-American-British 



classification M5) does not have a worse prognosis than other subtypes of 
Acute Myeloid Leukemia: a report from the eastern cooperative oncology 
group. J Clin Oncol 2004, 22:1 276-1 286. 

28. Bullinger L, Dohner K, Bair E, Frohling S, Schlenk RF, Tibshirani R, Dohner H, 
Pollack JR: Use of gene-expression profiling to identify prognostic 
subclasses in adult acute myeloid leukemia. N Engl J Med 2004, 
350:1605-1616. 

29. Schoch C, Schnittger S, Klaus M, Kern W, Hiddemann W, Haferlach T: AML 
with 1 1q23/M/l abnormalities as defined by the WHO classification: 
incidence, partner chromosomes, FAB subtype, age distribution, and 
prognostic impact in an unselected series of 1897 cytogenetically 
analyzed AML cases. Blood 2003, 102:2395-2402. 

30. Parker H, Rose-Zerilli MJJ, Parker A, Chaplin T, Wade R, Gardiner A, Griffiths M, 
Collins A, Young BD, Oscier DG, Strefford JC: 13q deletion anatomy and 
disease progression in patients with chronic lymphocytic leukemia. 
Leukemia 201 1,25:489-497. 

31. Bibikova M, Barnes BT, Chan T, Ho V, Klotzle B, Le JM, Delano D, Zhang L, 
Schroth GP, Gunderson KL, Fan J-B, Shen R: High density DNA methylation 
array with single CpG site resolution. Genomics 201 1, 98:288-295. 

32. Sarova I, Brezinova J, Zemanova Z, Bystricka D, Krejcik Z, Soukup P, Vydra J, 
Cermak J, Jonasova A, Michalova K: Characterization of chromosome 1 1 
breakpoints and the areas of deletion and amplification in patients with 
newly diagnosed acute myeloid leukemia. Genes Chromosomes Cancer 
2013,52:619-635. 

33. Calin GA, Cimmino A, Fabbri M, Ferracin M, Wojcik SE, Shimizu M, Taccioli C, 
Zanesi N, Garzon R, Aqeilan Rl, Alder H, Volinia S, Rassenti L, Liu X, Liu C-G, 
Kipps TJ, Negrini M, Croce CM: MiR-15a and miR-16-1 cluster functions in 
human leukemia. Proc Natl Acad Sci U S A 2008, 105:5166-5171. 

34. Morenos L, Saffery R, Mechinaud F, Ashley D, Elwood N, Craig JM, Wong 
NC: Evaluation of MicroRNA expression in patient bone marrow aspirate 
slides. PLoS One 2012, 7:e42951. 

35. Monteys AM, Spengler RM, Wan J, Tecedor L, Lennox KA, Xing Y, Davidson 
BL: Structure and activity of putative intronic miRNA promoters. RNA 
2010, 16:495-505. 

36. Corcoran MM, Hammarsund M, Zhu C, Lerner M, Kapanadze B, Wilson B, 
Larsson C, Forsberg L, Ibbotson RE, Einhorn S, Oscier DG, Grander D, 
Sangfelt O: DLEU2 encodes an antisense RNA for the putative bicistronic 
RFP2/LEU5 gene in humans and mouse. Genes Chromosomes Cancer 2004, 
40:285-297. 

37. Liu Q, Fu H, Sun F, Zhang H, Tie Y, Zhu J, Xing R, Sun Z, Zheng X: miR-16 
family induces cell cycle arrest by regulating multiple cell cycle genes. 
Nucleic Acids Res 2008, 36:5391-5404. 

38. Aqeilan Rl, Calin GA, Croce CM: miR-15a and miR-16-1 in cancer: discovery, 
function and future perspectives. Cell Death Differ 2010, 17:215-220. 

39. Cittelly DM, Das PM, Salvo VA, Fonseca JP, Burow ME, Jones FE: Oncogenic 
HER2D16 suppresses miR-15a/16 and deregulates BCL-2 to promote 
endocrine resistance of breast tumors. Carcinogenesis 2010, 31:2049-2057. 

40. Hanlon K, Rudin CE, Harries LW: Investigating the targets of MIR-15a and 
MIR-16-1 in patients with Chronic Lymphocytic Leukemia (CLL). PLoS One 
2009, 4:e7169. 

41. Roccaro AM, Sacco A, Thompson B, Leleu X, Azab AK, Azab F, Runnels J, 
Jia X, Ngo HT, Melhem MR, Lin CP, Ribatti D, Rollins BJ, Witzig TE, Anderson 
KC, Ghobrial IM: MicroRNAs 15a and 16 regulate tumor proliferation in 
multiple myeloma. Blood 2009, 113:6669-6680. 

42. Zhao H, Kalota A, Jin S, Gewirtz AM: The c-myb proto-oncogene and 
microRNA-15a comprise an active autoregulatory feedback loop in 
human hematopoietic cells. Blood 2009, 113:505-516. 

43. Reutershan J, Cagnina RE, Chang D, Linden J, Ley K: Therapeutic anti- 
inflammatory effects of myeloid cell adenosine receptor A2a stimulation in 
lipopolysaccharide-induced lung injury. J Immunol 2007, 1 79:1 254-1 263. 

44. Wolff F, Loipetzberger A, Gruber W, Esterbauer H, Aberger F, Frischauf AM: 
Imiquimod directly inhibits Hedgehog signalling by stimulating 
adenosine receptor/protein kinase A-mediated GLI phosphorylation. 
Oncogene 2013, 32:5574-5581. 

45. Andersson A, Ritz C, Lindgren D, Eden P, Lassen C, Heldrup J, Olofsson T, Rade J, 
Fontes M, Porwit-MacDonald A, Behrendtz M, Hoglund M, Johansson B, Fioretos 
T: Microarray-based classification of a consecutive series of 121 childhood 
acute leukemias: prediction of leukemic and genetic subtype as well as of 
minimal residual disease status. Leukemia 2007, 21:1 198-1203. 

46. Neff T, Armstrong SA: Recent progress toward epigenetic therapies: 
the example of mixed lineage leukemia. Blood 2013, 121:4847-4853. 



Morenos et al. Molecular Cancer 2014, 13:123 
http://www.molecular-cancer.eom/content/1 3/1/1 23 



Page 1 1 of 1 1 



47. 



49. 



50. 



51. 



52. 



53. 



54. 



55. 



56. 



57. 



58. 



59. 



60. 



62. 



63. 



64. 



Haferlach T, Schoch C, Schnittger S, Kern W, Loffler H, Hiddemann W: 
Distinct genetic patterns can be identified in acute monoblastic and 
acute monocytic leukaemia (FAB AML M5a and M5b): a study of 124 
patients. Br J Haematol 2002, 1 1 8:426-431 . 

Meyer C, Kowarz E, Hofmann J, Renneville A, Zuna J, Trka J, Abdelali RB, 
Macintyre E, De Braekeleer E, De Braekeleer M, Delabesse E, de Oliveira MP, 
Cave H, Clappier E, van Dongen JJM, Balgobind BV, van den Heuvel-Eibrink MM, 
Beverloo HB, Panzer-Grumayer R, Teigler-Schlegel A, Harbott J, Kjeldsen E, 
Schnittger S, Koehl U, Gruhn B, Heidenreich 0, Chan LC, Yip SF, Krzywinski M, 
Eckert C, et al: New insights to the MLL recombinome of acute leukemias. 
Leukemia 2009, 23:1490-1499. 

Schoch C, Kohlmann A, Dugas M, Kern W, Schnittger S, Haferlach T: Impact 
of trisomy 8 on expression of genes located on chromosome 8 in 
different AML Subgroups. Genes Chromosomes Cancer 2006, 45:1 164-1 168. 
Wong NC, Ashley D, Chatterton Z, Parkinson-Bates M, Ng HK, Halemba MS, 
Kowalczyk A, Bedo J, Wang Q, Bell K, Algar E, Craig JM, Saffery R: A distinct 
DNA methylation signature defines pediatric pre-B cell acute lymphoblastic 
leukemia. Epigenetics 2012, 7:535-541. 

Rosenfeld C, Goutner A, Choquet C, Venuat AM, Kayibanda B, Pico JL: 

Phenotypic characterization of a unique non-T, non-B acute lymphoblastic 

leukaemia cell line. Nature 1977, 267:841-843. 

Kaplan J, Mastrangelo R, Peterson WDJ: Childhood lymphoblastic 

lymphoma, a cancer of Thymus-derived lymphocytes. Cancer Res 1974, 

34:521-525. 

Asou H, Tashiro S, Hamamoto K, Otsuji A, Kita K, Kamada N: Establishment 

of a human acute myeloid leukemia cell line (Kasumi-1) with 8;21 

chromosome translocation. Blood 1991, 77:2031-2036. 

Tsuchiya S, Yamabe M, Yamaguchi Y, Kobayashi Y, Konno T, Tada K: 

Establishment and characterization of a human acute monocytic 

leukemia cell line (THP-1). Int J Cancer 1980, 26:171-176. 

Lange B, Valtieri M, Santoli D, Caracciolo D, Mavilio F, Gemperlein I, Griffin C, 

Emanuel B, Finan J, Nowell P: Growth factor requirements of childhood 

acute leukemia: establishment of GM-CSF-dependent cell lines. Blood 

1987, 70:192-199. 

Du P, Kibbe WA, Lin SM: lumi: a pipeline for processing lllumina 
microarray. Bioinformatics 2008, 24:1547-1548. 

Smyth GK: limma: Linear Models for Microarray Data. New York, USA: 
Springer; 2005. 

Eden E, Navon R, Steinfeld I, Lipson D, Yakhini Z: GOrilla: a tool for 
discovery and visualization of enriched GO terms in ranked gene lists. 

BMC Bioinformatics 2009, 10:48. 

Li J, Smyth P, Flavin R, Cahill S, Denning K, Aherne S, Guenther SM, O'Leary 
JJ, Sheils O: Comparison of miRNA expression patterns using total RNA 
extracted from matched samples of formalin-fixed paraffin-embedded 
(FFPE) cells and snap frozen cells. BMC Biotechnol 2007, 7:36. 
Borze I, Guled M, Musse S, Raunio A, Elonen E, Saarinen-Pihkala U, 
Karjalainen-Lindsberg M-L, Lahti L, Knuutila S: MicroRNA microarrays on 
archive bone marrow core biopsies of leukemias — Method validation. 
Leuk Res 2011,35:188-195. 

Livak KJ, Schmittgen TD: Analysis of relative gene expression data using 
real-time quantitative PCR and the 2" DDCt method. Methods 2001, 25:402-408. 
Betel D, Wilson M, Gadow A, Merks DS, Sander C: The microRNA.org 
resource: targets and expression. Nucleic Acids Res 2008, 36:D149-D153. 
John B, Enright AJ, Aravin A, ETuschl T, Sander C, Merks DS: Human 
microRNA targets. PLoS Biol 2004, 2:e363. 

Dweep H, Sticht C, Pandey P, Gretz N: miRWalk-database: prediction of 
possible miRNA binding sites by "walking" the genes of three genomes. 

J Biomed Inform 201 1, 44:839-847. 



doi:1 0.1 1 86/1 476-4598-1 3-1 23 

Cite this article as: Morenos et al.: Hypermethylation and down-regulation 
of DLEU2 in paediatric acute myeloid leukaemia independent of embedded 
tumour suppressor miR-15a/16-1. Molecular Cancer 2014 13:1 23. 



Submit your next manuscript to BioMed Central 
and take full advantage of: 

• Convenient online submission 

• Thorough peer review 

• No space constraints or color figure charges 

• Immediate publication on acceptance 

• Inclusion in PubMed, CAS, Scopus and Google Scholar 

• Research which is freely available for redistribution 



Submit your manuscript at 
www.biomedcentral.com/submit 



o 



BioMed Central 



